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THE  EFFECT  OF  BEAM  QUALITY  ON  THE 
FREE-ELECTRON  LASER 


I  I.  INTRODUCTION 


One  of  the  most  important  current  problems  involved  in  the  design  and  operation  of  the 
Free-Electron  Laser  (FEL)  is  the  effect  of  an  initial  thermal  spread  in  the  electron  beam  ve¬ 
locity  distribution  upon  both  the  linear  gain  and  nonlinear  saturation  efficiency  of  the  interac- 
}  tion.  The  physical  mechanism  underlying  the  concept  of  the  FEL  is  the  axial  bunching  of  an 

I  electron  beam  in  the  ponderomotive  wave  formed  by  the  beating  of  radiation  and  wiggler 

'•  fields,  and  maximum  operational  efficiencies  are  obtained  when  saturation  occurs  by  the  trap- 

I  ping  of  the  bulk  of  the  electron  beam  in  the  ponderomotive  potential.  Such  a  process 

s,  requires  a  highly  coherent  interaction  involving  both  a  narrow  radiation  spectrum  and  an 

j  extremely  low  axial  velocity  spread  in  the  electron  beam.  Observe  that  a  monoenergetic  elec- 

L  tron  distribution  is  not  sufficient  to  ensure  operation  in  the  trapped-particle  regime,  and  the 

I  partition  of  energy  between  the  axial  and  transverse  directions  is  important. 


The  motivation  of  the  present  work  is  to  develop  a  fully  three-dimensional  nonlinear 
formulation  of  the  FEL  amplifier  for  a  configuration  which  consists  of  an  electron  beam  pro¬ 
pagating  through  a  loss-free  cylindrical  waveguide  in  the  presence  of  a  helically  symmetric 
wiggler  and  uniform  axial  guide  magnetic  field.  The  analysis  we  employ  is  an  extension  of 
previously  described  nonlinear  analysis1-2  which  treated  the  problem  in  the  absence  of  an  ini¬ 
tial  beam  velocity  spread.  The  formulation  involves  the  derivation  and  numerical  solution  of 
a  set  of  coupled  nonlinear  differential  equations  which  describe  the  self-consistent  evolution 
of  both  the  trajectories  of  an  ensemble  of  electrons  and  the  electromagnetic  fields.  Space- 
charge  fields,  however,  are  neglected  and  the  analysis  is  restricted  to  the  high-gain  Compton 
(or  "strong-pump")  regime  of  operation.  The  adiabatic  injection  of  the  electron  beam  is 
modeled  by  allowing  the  wiggler  field  amplitude  to  increase  slowly  from  zero  to  a  constant 
level.  The  initial  electron  beam  distribution  (i.e.,  external  to  the  wiggler  field  region)  is 
assumed  to  be  monoenergetic  but  of  finite  emittance.3  Since  we  include  the  presence  of  an 
axial  guide  magnetic  field  in  the  formulation,  the  concept  of  beam  emittance  is  ambiguous 
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and  can  be  defined  in  various  ways.  For  convenience  in  this  paper,  however,  we  use  the 
term  merely  to  describe  an  electron  beam  which  is  characterized  by  an  initial  momentum 
spread  in  (pL ,  p. )  that  describes  a  ring  distribution  subject  to  the  condition  that  p‘  -  p?  = 
constant.  Such  a  distribution  describes  an  electron  beam  with  both  a  transverse  and  axial 
energy  spread,  and  the  magnitude  of  the  spread  provides  a  measure  of  the  electron  beam 
quality  for  FEL  applications.  The  propagation  of  the  beam  through  the  adiabatic  entry  region 
of  the  wiggler,  therefore,  describes  the  self-consistent  evolution  of  the  initial  distribution  due 
to  the  effect  of  the  wiggler  field.  It  should  be  remarked  that  there  is  no  fundamental  diffi¬ 
culty  in  the  inclusion  of  a  spread  in  the  total  energy  as  well,  and  the  formulation  can  be 
readily  generalized  to  describe  a  spread  in  the  total  energy.  This  occurs,  however,  at  the 
expense  of  the  inclusion  of  a  larger  number  of  electrons  in  the  numerical  simulation. 

The  organization  of  the  paper  is  as  follows.  The  general  equations  are  presented  in  Sec. 
II.  The  coupled  particle  and  field  equations  are  solved  numerically,  and  the  results  are 
described  in  Sec.  Ill  for  a  variety  of  parameters  in  order  to  illustrate  the  effect  of  the  momen¬ 
tum  spread  on  the  operational  efficiency  of  the  FEL  amplifier.  Finally,  a  summary  and  con¬ 
clusion  is  given  in  Sec.  V. 


II.  THE  GENERAL  EQUATIONS 


The  physical  configuration  we  employ  describes  the  propagation  of  an  electron  beam 
through  a  loss-free  cylindrical  waveguide  in  the  presence  of  an  axial  guide  field  and  a  helical 
wiggler  field.  The  analysis  is  three-dimensional  and  we  represent  the  static  magnetic  fields  in 
cylindrical  coordinates  as 


Bj(xi  «=  2?oe,  +  2Bk(i) 


/{(k)e,  cos  \  — “  /j(X)  e#  sin  x  +  A  (A)  e.  sin  x 
\ 


(1) 


where  B0  and  B „  represent  the  axial  and  wiggler  field  amplitudes,  k  =  A„r,  *  =  0  -  kw r, 
kH  =  2i t/a„.  (where  \M  is  the  wiggler  period),  and  I„ (k)  and  (A >  are  the  modified  Bessel 
function  and  its  derivative.  In  order  to  model  the  realistic  injection  of  the  electron  beam,  we 
allow  the  wiggler  amplitude  to  increase  adiabatically  to  a  constant  amplitude  over  A*  wiggler 
periods,  and  write 


BH  sin2  (kwz/ 4A'm.)  ,  0 
B»  ,  r  >  A',Am 


(2) 


Space-charge  fields  are  neglected,  and  the  boundary  conditions  of  the  radiation  field  at  the 
waveguide  wall  are  satisfied  by  expanding  the  vector  potential  in  terms  of  the  basis  functions 
of  the  vacuum  guide.  We  restrict  the  analysis  to  the  TE  modes  and  write 
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8A (x.t)  -  — — J;(kinr)  er  sin  a,  +  J{ (k,„r)  e„  cos  a,  ,  (3) 

,1—  l 

where  for  frequency  <o  and  wavenumber  JtG) 

a /  =  J*0  dz'k(z')  +  19  —  (ut ,  (4) 

£/„  s  /Rf  (  Rg  =  waveguide  radius)  for  J't(x'!n)  -  0,  and  /,Cr)  and  7',.  (x)  are  the  regular 
Bessel  function  of  the  first  kind  and  its  derivative.  Both  the  mode  amplitude  8A,„(z)  and  the 
wavevector  k(z)  are  assumed  to  vary  slowly  over  a  wave  period.  The  electron  beam  is 
assumed  to  be  initially  (i.e.,  at  z  -  0)  monoenergetic,  but  with  a  momentum  distribution  in 
the  form  of  a  "ring”  in  the  initial  p10  and  p.0  as  follows: 

foW  "  A  exp  (-(/?: o  -  p0)2/\p?)  8  (pi  -  Plo  ~  P:2o)Q(p-.o)  (5) 

where  (p0.  Ap.)  represent  the  total  momentum  and  axial  momentum  spread  of  the  beam, 
Q(x )  is  the  Heaviside  function,  and  the  normalization  constant  is 


A  -  (it  J*/%r0  exp  (-  (p. o  -  po^/Apz))  . 


Note,  therefore,  that  Ip.  describes  a  spread  in  both  axial  and  transverse  momenta. 

The  derivation  of  the  equations  governing  the  evolution  of  8A,„(z )  and  k(z)  proceeds  in 
the  same  manner  described  by  Ganguly  and  Freund.1-2  and  we  find  that 


d2  u2  V,  7T’  +  V,  »y*> 

—  8a,n+  *j-k--krH  8aln  —  -j-  0.o  Ht„  <J— ,  ~  >■ 


I 


2 kU2  ±(k"28aj  -  ~  0zO  H,n  < 


V[  -  v,  r/-1 


where  8a,„  s  e8/J;„  /'me2,  /3;0  =  v;0/c.  w42  s  4m e2nb/m ,  (v,,  v;)  are  the  transverse  components 
of  the  velocity  relative  to  the  basis  vectors  e,  -  eT  cos  kvz  +  ev  sin  /twr, 

—  -eT  sin  kvz  -  e„  cos  /c„z. 


^  (■*'£  -  /2)  7,20r'J 


r/->  =  F/-’  sin  i/r,  +  (7/=l  cos  <i»,,  andH,'/=>  =  cos  \p,  -  GilS)  sin  > h, .  In  addition, 
F/=J  S  7;_ \(ki„r)  cos  (/  -  1)  x  —  cos  <7+l)  x. 


=  < Ar,w r *  sin  f /—  l )  <  =  '<,„/•>  sin  '/-I/  x- 
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is  the  particie  phase  relative  to  the  ponderomotive  wave  -  mt0).  These  equations  are 
formally  identical  to  those  derived  previously,  except  that  average  operator  ( Ag  -  irRf) 


<("")>  =  \~2  L  d(t> 0  fo°  dP:0  T2-  exp  (-  (p.  o  -  Pq)2/±p{) 


4 irzRgl  Jo 


P  o 


x  f_v  dtb<\  <r ;,((i»0)  f  f  dffqdrqfq  cr_  (/-0,  0O)  (•••■) 


(12) 


now  includes  the  initial  momentum  space  coordinates  p:0  and  <t>0  =  ian~l(py0  /px0).  Note  that 
a-  <0O)  and  <r_{r0,  Qq)  describe  the  initial  electron  distribution  in  phase  and  cross-section.  For 
simplicity,  we  assume  the  beam  to  be  uniformly  distributed  in  phase  (cr,  «■  1)  and  cross- 
section  for  a  given  beam  radius  Rb  (i.e.,  <rL  -  1  for  0  <  r  <  Rb). 


Finally,  in  order  to  complete  the  formulation,  the  electron  orbit  equations  must  be  speci¬ 
fied  by  the  Lorentz  force  equations.  For  the  particular  static  and  fluctuating  fields  described 
previously,  we  find  that  for  the  TE,„  mode 


v.~pi  -  — -[fto-  +  2flw/,(x)  sinxlp-*  +  —  12„p./2(x)  sin  2x 

'  dz  y  '  y 


— \-mc5a,nl(<D  ~  Arv.)  W{~'  -  2k, n w2J,  (k,„r)  cos  a, 

4m 

-r  tov.m, 

v.  —rP"  -  —  [Slq  -  v£„v.  J-  212, Ii(\)  sin  xl/?i 
■  az  •/ 


(13) 


- -!),;.[/o(l)  -*■  /;  (X )  COS  2x1 

y 


t  —rnc&a,n[(<n  —  kv.)Tt'~'  -  Ik^v^JAk^r)  cos  at 

+  r,>v.  wr'\. 


(14) 


v.-7-p.  -  — n,p-(/o(X)  4-  /-.(X)  cos  2x1  — -12 ,^i/;(X)  sin  2x 
■«•-/•  •  y 


-  —mc8a,„  (fc(v,  tF/-*  -  v:r/_l) 


r (v:  r/*'  + 


(15) 
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where  f20.»  —  leB^Jmc i,  y  =  (1  -  v;/c:)-!/\  and  r.„  =  d(tn8a,n)/dz.  In  addition. 


v.-f-.r  —  v,  cos  kw:  —  v,  sin  kv r, 
'  az  ' 


1,1 1 

$ 

$ 

♦%* 

w 


fft I 

If 

:|a 

tjy 

s 


a 


v.  — y  —  Vj  sin  cos  £v  j, 
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III.  NUMERICAL  SIMULATION 

We  first  consider  several  cases  corresponding  to  the  low  axial  guide  field  regime  for 
which  fl0  <  vo^v„.  Case  1  is  that  of  a  35-GHz  amplifier,  and  we  assume  that  B0-  1-3  kG, 
Bw  -  2  kG,  \w  -  1.175  cm,  and  ;V„  -  10.  The  electron  beam  is  characterized  by  an  energy  of 
250  keV,  a  current  of  35  A,  and  an  initial  radius  Rb  -  0.155  cm.  This  case  has  been  con¬ 
sidered  previously  in  the  absence  of  an  initial  velocity  spread1,  and  a  maximum  efficiency 
tj0  -  21.4%  was  found  for  a  frequency  at /ckw  —  1.3  (33.2  GHz)  in  the  TEn  mode  (Rg  —  0.36626 
cm).  The  effect  of  an  initial  velocity  spread  is  shown  in  Fig.  1,  in  which  we  plot  the  satura¬ 
tion  efficiency  versus  lp./p0.  The  decrease  in  efficiency  with  ip.  is  quite  rapid,  and  an  order 
of  magnitude  decrease  in  tj  is  found  to  occur  as  the  velocity  spread  increases  to  lp./p0  —  2%. 
In  terms  of  an  axial  energy  spread  we  write 

1-  i  +  lysBo2—  .  (19) 

yo  Po 

where  y0  -  (1  +  p$/m:cz)l/2,  and  /30  -  po/y^mc.  Hence,  this  order  of  magnitude  decrease  in 
the  efficiency  corresponds  to  Sy:/yo  —  2.6%. 

As  shown  by  the  numerical  results,  the  saturation  efficiency  is  highly  sensitive  to  the 
axial  energy  spread  of  the  beam.  Typically,  the  axial  energy  spread  must  be  much  less  than 
the  depth  of  the  ponderomotive  potential  in  order  for  the  trapping  mechanism  to  be  effective. 
On  the  basis  of  perturbations  about  an  idealized  three-dimensional  class  of  helical  trajec¬ 
tories,4  we  obtain  an  estimate  of  the  effective  ponderomotive  "depth"  of 


in .  4  /sE*Jl ta,  (2o> 
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where  (v„,  v.,)  are  the  transverse  and  axial  velocity  components  of  the  helical  orbits, 
v  -  -  1 1  —  v -  /c-r‘,  r0  is  the  radius  of  the  orbit  about  the  axis  of  symmetry.  j>  :s  a  complicated 


v '-'A 


function  of  the  orbit  parameters  given  in  Ref.  2,  and  the  normalized  radiation  amplitude  Ha,„ 
is  to  be  evaluated  at  saturation  (in  the  absence  of  thermal  spread).  Eq.  (20)  provides  only  a 
crude  estimate  of  the  trapping  depth  since  the  trajectories  both  in  simulation  and  in  actual 
experiments  can  differ  from  the  heiical  orbits.  However,  if  we  relate  Avl;/v-  —  A p./p0  and  use 
Eq.  (20)  to  obtain  a  measure  of  the  effective  ponderomotive  energy  spread  Aypon<j/'-y0.  then  we 
find  that  Ay^  y0  —  21%  for  parameters  associated  with  Case  1.  This  indicates  that  the  effi¬ 
ciency  has  decreased  by  an  order  of  magnitude  as  the  axial  energy  spread  approaches  about 
10%  of  the  ponderomotive  "trapping  depth." 

Case  2  corresponds  to  an  FEL  experiment  at  the  Naval  Research  Laboratory5-0  which 
makes  use  of  an  electron  beam  with  an  energy  of  750  keV,  a  current  of  200  A.  and  an  initial 
beam  radius  Rb  *  0.5  cm.  The  experiment  employs  a  helical  wiggler  with  Bw  —  1  kG,  * 
4  cm,  and  .V„  —  6.  No  axial  guide  field  is  imposed.  This  case  has  also  been  studies  previ¬ 
ously,2  and  an  efficiency  of  170  =  7.27%  was  found  at  a  frequency  cu /ckw  =  8.3  (62.3  GHz)  in 
the  TE\\  mode  (Rg  -  1.5  cm).  The  saturation  efficiency  is  shown  in  Fig.  2  as  a  function  of 
Ap./poi  and  it  is  evident  that  an  order  of  magnitude  decrease  in  the  efficiency  has  occurs  as 
the  momentum  spread  increases  to  A p./p0  — ;  0.21%.  This  corresponds  to  an  axial  energy 
spread  of  Ay./yo  —  1%.  For  this  case,  the  ponderomotive  trapping  depth  is  Ay^yi)  —  15% 
and.  as  for  Case  1,  the  efficiency  decreases  by  an  order  of  magnitude  as  the  axial  energy 
spread  approaches  10%  of  the  "trapping  depth."  Finally,  we  interpret  the  relatively  sharp  tran¬ 
sition  of  the  scaling  of  77  with  A p.  which  occurs  for  Ap./po  —  0.1-0.15  as  due  to  the  transition 
between  the  cold  and  thermal  beam  regimes.  This  will  be  discussed  in  more  detail  later  in 
this  section.  We  remark  that  such  rapid  declines  in  the  efficiency  with  increasing  momentum 
spread  as  shown  in  Cases  1  and  2  are  consistent  with  results  found  previously  by  Kwan  and 
Snell7  by  means  of  a  one-dimensional  particle- in-ceil  simulation  code. 

In  Case  3.  we  investigate  the  effect  of  a  relatively  small  axial  guide  field  on  the  scaling  of 
the  efficiency  for  parameters  consistent  with  Case  2  by  inclusion  of  a  guide  field  of  Bn  -  1  kG. 
This  is  still  far  below  the  resonance  condition  for  n„  —  y^v,,;  however,  the  axial  velocity  is 
sufficiently  reduced  that  the  frequency  regime  over  which  gain  occurs  drops  and  we  find  an 
efficiency  in  the  absence  of  a  momentum  spread  770  —  8.37%  at  a  frequency  (aidcw  =  7.5 
(56.2  GHz).  The  scaling  of  the  efficiency  versus  Ap./p0  for  this  case  is  shown  in  Fig.  3.  It  is 
clear  from  the  figure  that  the  effect  of  the  guide  field  has  been  to  reduce  the  sensitivity  of  the 
efficiency  to  the  momentum  spread  since  the  efficiency  has  decreased  by  only  about  80%  rela¬ 
tive  to  the  zero  momentum  spread  result  for  Ap:/p0  =  0.6%.  This  contrasts  with  Case  2  in 
which  an  order  of  magnitude  decrease  in  tj  was  observed  for  Ap./po  —  0.21%.  We  also  note 


that  lp./po  —  0.6%  corresponds  to  an  axial  energy  spread  of  ly./y 0  =  2.9%.  Since  Av?ona,v0 
—  19.1%  for  this  case,  the  80%  decrease  in  r\  corresponds  to  an  axial  energy  spread  of  about 
15%  of  the  "trapping  depth." 

We  now  address  the  effect  of  a  strong  axial  field  on  the  sensitivity  of  the  efficiency  to 
the  beam  emittance.  We  observe  that  the  effective  ponderomotive  energy  spread  (20)  will 
increase  as  30  approaches  the  resonance  at  fl 0  —  ?o^h.v!|.  because  vw,  <4> ,  and  5  are 
enhanced  near  the  resonance  while  vt|  declines.  Operation  close  to  resonance,  therefore, 
should  result  in  a  decreased  sensitivity  to  beam  emittance.  In  order  to  verify  this  conclusion 
numerically,  we  consider  the  effect  of  finite  emittance  in  the  regime  in  which  n0  ^  yo^vM- 
The  example  (Case  4)  of  interest  corresponds  to  the  propagation  of  a  1  MeV,  50  A  beam 
with  Rb  -  0.2  cm  through  a  drift  tube  of  radius  0.5  cm.  The  static  magnetic  fields  are: 
30  -  U.15kG,  Bw  -  l.OkG,  \w  -  3.0cm,  and  Nw  -  10.  This  case  exhibits  a  high  efficiency  and 
broad  bandwidth  in  the  zero  momentum  spread  limit  with  a  peak  efficiency  of  rj0  =46.9%  at 
eo/c  kv  -  3.0  (30  GHz)  in  the  TEU  mode.1  The  effect  of  finite  emittance  is  shown  in  Fig.  4  in 
which  we  plot  the  normalized  emittance  versus  IpJpo-  Clearly,  the  efficiency  decreases  with 
emittance  far  more  slowly  than  in  either  of  the  other  cases,  and  decreases  by  an  order  of 
magnitude  for  \p./p0  =  5-6%  (corresponding  to  \y.  /y0  ~  25%).  As  a  result,  we  conclude 
that  dramatic  decreases  in  the  sensitivity  of  the  efficiency  upon  the  initial  beam  emittance  can 
be  obtained  from  operation  with  strong  axial  fields  close  to  the  resonance  at  n0  3*  y,  kv  vu. 

Although  the  effect  is  more  evident  in  Cases  2-4,  all  of  the  cases  shown  exhibit  a  rela¬ 
tively  sudden  change  in  the  scaling  of  the  efficiency  versus  \p:/p0  which  we  attribute  to  the 
transition  from  the  cold  beam  limit  to  the  "thermal”  regime  in  which  the  wave  is  resonant 
with  only  a  small  fraction  of  the  electron  beam  due  to  finite  momentum  spread  effects.  An 
estimate  of  where  this  transition  occurs  may  be  obtained  by  noting  that  the  cold  beam  limit 
occurs  for  \lmk!  Rek\  »  lpz/po-  We  note  that  the  ratio  \Imk/Rek\  -  3.5%,  0.48%,  0.62%,  and 
1.8%  for  Cases  1-4  respectively.  The  determination  of  the  precise  range  of  the  transition 
range  between  the  cold  and  thermal  regimes  requires  a  detailed  analysis  of  the  linear  disper¬ 
sion  properties  of  the  modes  in  three-dimensions.  However,  we  expect  that  this  transition 
occurs  for  \p./p0  <  \Imk/ Rek\  <  3Ap./p0,  which  provides  reasonable  agreement  with  the 
numerical  results. 

IV.  SUMMARY  AND  CONCLUSION 

In  this  paper  we  have  developed  a  three-dimensional  nonlinear  theory  and  simulation  of 
the  FEL  amplifier  in  the  high-gain  Compton  regime  with  the  inclusion  of  a  finite  electron 


beam  momentum  spread.  The  formulation  describes  the  adiabatic  injection  of  the  eiectron- 
beam  into  the  wiggler  field  region;  hence  the  evolution  of  the  beam  distribution  'initially 
defined  external  to  the  wiggler)  as  the  beam  enters  the  wiggler  is  described  in  a  self- 
consistent  manner.  Results  of  the  simulation  indicate  that  in  the  low  axial  guide  field  regime 
the  initial  axial  energy  spread  must  be  less  than  approximately  10%  of  the  ponderomotive 
trapping  depth  at  saturation  (in  the  absence  of  momentum  spread)  in  order  for  high  efficiency 
operation  to  be  achieved.  However,  the  effect  of  the  axial  guide  field  appears  to  decrease  the 
sensitivity  of  the  efficiency  to  the  momentum  spread.  This  occurs  due  to  the  resonant 
enhancement  of  the  ponderomotive  trapping  depth  (20),  and  is  most  pronounced  for  high 
axial  guide  fields  near  the  resonance  at  n0  <  yt)kwv\\.  Finally,  a  relatively  sharp  transition  in 
the  scaling  of  the  efficiency  with  the  axial  momentum  spread  is  observed  in  simulation  which 
is  attributed  to  the  transition  from  the  sold  beam  to  the  "thermal"  beam  regime. 
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Fig.  2  —  Graph  of  the  efficiency  versus  axial  momentum  spread  for  Case  2 
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